Viscoelastic properties of attractive and repulsive colloidal glasses 
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We report a numerical study of the shear viscosity and the frequency dependent elastic moduli 
close to dynamical arrest for a model of short-range attractive colloids, both for the repulsive and 
the attractive glass transition. Calculating the stress autocorrelation functions, we find that density 
fluctuations of wavevectors close to the first peak in the structure factor control the viscosity rise on 
approaching the repulsive glass, while fluctuations of larger wavevectors control the viscosity close 
to the attractive glass. On approaching the glass transition, the viscosity diverges with a power law 
with the same exponent as the density autocorrelation time. 

PACS numbers: 82.70.Dd, 82.70.Gg, 64.70.Pf 
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Colloidal dispersions, in addition to their technological 
relevance, play an important role in the development of 
basic physical sciences, in particular in the fascinating 
field of formation of disordered arrested states, glasses 
and gels. The possibility of tailoring shape, size and 
structure of the colloidal particles makes it possible to 
design specific colloidal interaction potentials and to use 
light scattering to directly probe even single particle mo- 
tion. Recently, the application of experimental[l|, |2j, |3( 
and theoretical toolsQ OJ ||J to the analysis of the glass 
transition in short-ranged attractive colloids, has shown 
an extremely rich scenario, with no analogue in atomic 
systems 0- The standard packing-driven hard-sphere 
glass transition transforms - discontinuously in some 
cases - into a novel-type of glass transition driven by the 
short-range attractions. The connection between this at- 
tractive transition and gelation is still matter of debate. 

Previous numerical work on short-range attractive col- 
loids has mostly focused on the behavior of single parti- 
cle diffusion and of the time dependence of density au- 
tocorrelation function 3> 9 (t) close to dynamical arrest. 
Despite the strong link with experiments and the rele- 
vance to industrial applications, the numerical evalua- 
tion of the viscosity, r\, and viscoelastic properties f)(uj) 
has not been performed, since significant computational 
effort is required for accurate calculation of r)(u:), even 
more for states close to dynamical arrest. The behavior 
of r\ close to dynamical arrest in colloidal systems is under 
theoretical and experimental investigation^, 0, 0, Hlf . 
Measurements in dense hard sphere colloids show a di- 
vergence of rj in the vicinity of the repulsive glass point, 
its exact location and functional form being still under 
debate 0,^3- For colloidal gels, a power law divergence 
has been reported in connection to the gel transition 
[l3|. The Mode Coupling Theory for glass transitions [li) 
(MCT) , the theoretical framework which anticipated the 
presence of novel dynamic phenomena in short-range at- 
tractive colloidsjj, M M predicts asymptotic power law 
divergence of rj with the distance to the transition, with 



identical exponents to the divergence of the time scale, 
r, and the inverse of the self diffusion coefficient 1/Do. 

In this Letter we report an extensive numerical study of 
the viscoelastic behavior of a colloidal system, approach- 
ing both the repulsive and the attractive glass transition. 
We calculate the stress autocorrelation function, as well 
as the viscosity and the elastic modulii. We find that, 
close to both transitions, rj diverges with a power law, 
with the same exponent as the density relaxation time, 
but different from that of 1/Dq. Moreover, we provide 
evidence that the rise of ij is controlled by density fluctua- 
tions of wavevectors around the first peak in the structure 
factor close to the repulsive glass, and by fluctuations of 
larger wavevectors close to the attractive glass. 

Our system is composed of 1000 particles interacting 
via a steep repulsive (r~ 36 ) potential complemented by 
the Asakura-Oosawa (AO) short-range attractive poten- 
tial The attraction strength is measured in units of 
the polymer volume fraction, <p p and density is reported 
as colloid volume fraction, <fi c = 4/37ra[a 2 + 5 2 ]p, with p 
the number density, and the range of the interaction is 
£ = 0.1a 0|- Dynamics in colloidal systems can be con- 
veniently described by Brownian dynamics (BD), which 
accounts for the presence of solvent (not treating the full 
hydrodynamic interactions). Unfortunately, BD simula- 
tions are effectively less efficient than Newtonian dynam- 
ics (ND) due to the slower microscopic time scale. Since 
as far as slow dynamics is concerned, BD and ND are 
equivalent |17Lll9| we mostly work with the more efficient 
ND but perform additional simulations with BD at high 
friction coefficient, 7 = 50 |20j . to confirm the indepen- 
dence on the short-time dynamics. Dynamics in the cho- 
sen model system have been studied previously [Isl llflj. 
In the absence of polymers (</> p = 0) and with no repul- 
sive barrier, on increasing the particle packing fraction 
a repulsive glass transition is observed at <f>~? ~ 0.594. 
At fixed (j> c — 0.40, on increasing <p p , an attractive glass 
transition is observed at </£ ~ 0.4265. MCT analysis 
for <3> 9 (t) and Do has been performed previously for both 
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FIG. 1: Stress correlation functions for ND (grey lines) and 
BD (black lines), time rescaled, for 4> p = 0.30 and 0.41 as la- 
beled. For both states, Cj 3 / 3 ' (t) has been rescaled onto C^J 3 (t) 
by a factor 6.25. Inset: Full scale stress correlation functions 
for ND and BD (without time rescaling). 
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The shear viscosity rj is given by the Green-Kubo relation 



dtC aa (t) 



dt 



a</3 



(1) 

which expresses rj as the integral of the correlation func- 
tion of the non-diagonal terms of the microscopic stress 
tensor, a af3 = Yh=\ mv iaVi/3 ~ r ' j ^ ijf3 V (r. 



ij h 

where V is the volume of the simulation box, Vi a is the 
a-th component of the velocity of particle i, and V' is 
the derivative of the total potential. 

The inset in Fig. ^shows C aa (t) for two state points 
with 4> c = 0.40: <p p = 0.30 and <p p = 0.41, both for ND 
and BD. In the case of ND, clear oscillations in C aa (t) 
at small time (t < 0.5) are observed, caused by motion 
in the attractive well. These oscillations are completely 
damped in the BD case, confirming that this part of 
the decay of the correlation function has a microscopic 
(dynamic-dependent) origin and that short time oscilla- 
tions are not to be expected in real colloidal systems. 
When 4> p is increased, an additional slow relaxation pro- 
cess is present. Interestingly, the BD C aa (t) can be 
time-rescaled to collapse at long t on ND data (main 
figure), confirming that the long time dynamics of the 
system does not depend on the microscopic dynamics, in 
agreement with previous studies in glass forming atomic 
systems 01 . Therefore, we use ND hereafter to analyse 
the long time decay of C aa (t) and the viscosity close to 
both transitions. 

Fig. |21 shows the stress correlation function, C aa (t), 
for different <f> c or <j) p values on approaching the repulsive 
(upper panel) and attractive (lower panel) glass transi- 
tions. Correlation functions have been superimposed by 
a suitable scaling of time, t„ (insets of Fig. On ap- 
proaching the transition, C CTCT (t) slows down by several 
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FIG. 2: Stress correlation functions for different states ap- 
proaching the glass transitions. Upper panel: Repulsive 
glass transition. From right to left: <p c — 0.55, 0.57, 0.58. 
Lower panel: Attractive glass transition. From right to left: 
<t> p = 0.405,0.41,0.415,0.42 and cj> c = 0.40. Dashed lines are 
& q (t) 2 for qa — 4 at <j> c = 0.58 (repulsive glass) and qa = 11 
at <f) p = 0.415 (attractive glass). Insets: Time scale from 
the stress correlator vs. time scale from density correlator at 
the wavevectors given above; Caaija) = 1 (repulsive glass), 
Caa(ja) = 5 (attractive glass) and 4> 9 (t,j) = f q /e. 



orders of magnitude, with no appreciable change in its 
shape at long-time. A stretched exponential can be fit- 
ted to the decay with an exponent /3 = 0.48 and /3 = 0.25 
for the repulsive and attractive glasses, respectively. As 
can be seen from EqQJ the slowing down of C aa {i) is 
responsible of the rise of r\. 

The scaling property of C aa (t) at long times is remi- 
niscent of the scaling properties of & q (t). Indeed, within 
MCT, C nrT (t) can be related to &q(t), by means of 

Mm. 



c aa (t) 



k B T 
60tt 2 



dqq 4 
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-] 2 



dq 



(2) 



where S q is the static structure factor. Accordingly, the 
4> c or <f) p dependence of r CT should be identical to the <f> c 
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or 4> p dependence of r q , the characteristic time of $ 9 : 
^qijq) = /<j/ e ; with f q the non-ergodicity parameter. 
The two insets of Fig. |21 show a plot of T a vs. r q , para- 
metric in <j) c or (pp. For both types of glass transitions, the 
time scales obtained by rescaling C aa (t) are proportional 
to that found rescaling & q (t) (continuous lines). 

The dominant contribution to the g-integral in Eq. |2 
for times in the a-decay of $ q (t), can be revealed study- 
ing the function q 4 [d In S q /dq f q ] . This function oscil- 
lates with S q (out of phase ir/2), but its envelope shows 
a maximum close to the nearest neighbour peak of S q for 
the repulsive glass transition (qa w 3.75), and at higher q 
for the attractive one (the exact value depending on the 
details of the interaction potential). For the AO poten- 
tial, using Percus-Yevick closure, the dominant wavevec- 
tor is qa m 13, although the distribution is very wide. 
Interestingly, the wavevectors driving both the repulsive 
and the attractive ideal MCT glass transitions, as ob- 
served from the density correlators, are also the dominant 
ones in the calculation of C CT(T (t). 

Guided by Eq|3 we evaluate $ g from the simulations 
and compare with C aa {t) for different wavevec- 

tors, allowing for a scaling factor only in the amplitude. 
The best agreement is obtained for qa = 4 in the repul- 
sive glass transition and qa — 11 in the attractive glass 
(dashed lines in Fig. [21 , in good agreement with MCT 
predictions. The comparison shows that the long time 
decay of the stress correlation function, within the accu- 
racy of our numerical results, is adequately described by 
the dominant wavevector in Eq. [21 23] . 

We now turn to the study of 77, which is more efficiently 
calculated numerically using the Einstein relation 24] : 



rj = 1- lim j(AA(tf 



(3) 



where AA(t) is the integral from to t of the three off- 
diagonal terms of the stress tensor. Fig. [3] shows the 
behavior of 77 approaching the repulsive and the attrac- 
tive glass transitions, as a function of the distance to the 
transition points, (/>f and 4>p [3113 For both the repul- 
sive and attractive glass transitions, 77 is well described 
by a power law. In the attractive case, 77 shows a min- 
imum for intermediate (j> p , arising from the competition 
of the two competing arrest mechanisms [H Ejf . 

Fig. also shows T q for qa — 3.5 (the nearest neighbour 
peak in S q ), Dq, C CT<7 (0) and tjDq . For both transitions, 
both T q and Dq are linear in log-log and hence can also be 
well described by power laws. While 77 ~ r qi as confirmed 
by the constant ratio 77/r g (not shown), the product 77D0 
shows a clear trend, indicating that the power-law ex- 
ponents are very similar for 77 and r q but different for 
Dq. The best- fit values are reported in Table I for both 
transitions. The difference between the exponents 7^ 
and 7d indicates that the Stokes-Einstein relationship 
breaks down not only in the form Dq — r, as already 
shown in many different systems, but also in the form 
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FIG. 3: Viscosity, 77, time scale at q — 3.5, diffusion coefficient 
(-Do x 10 3 ) and Co-<r(0) for the repulsive glass transition (upper 
panel) and the attractive one (lower panel), as labeled. Note 
that 77D0 diverges in the transition point. The lines are the 
power law fittings to rj, r q and Do (see Table I). 
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TABLE I: Exponents obtained in fitting 77, r and 1/Do to 
(</>5? — C ) -7 , for the repulsive glass and (<f>p — (f> p ) ' for the 
attractive one. (*) The last point, 4> p — 0.42, is not considered 
in the fitting, due to the large error in determining 77. 



Da — rj. This breakdown, much more evident in the at- 
tractive glass case, is not consistent with MCT, although 
the theory correctly predicts the similarity between the 
exponents of T q and 77. The power-law divergence of 77 at 
the attractive glass transition agrees with recent experi- 
mental observations in colloidal gels ^ is also con- 
sistent with recent interpretation of data for the repulsive 
glass transition 0] > although this analysis is not conclu- 
sive, since an exponential divergence around <j> m = 0.64 
has also been suggested ^(J- The constant C CT(T (0) is the 
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FIG. 4: Elastic (upper panel) and viscous (lower 
panel) modulus for different states approaching the at- 
tractive glass transition: From right to left <f) p = 
0.30, 0.35, 0.375, 0.39, 0.40, 0.405, 0.41, 0.415, 0.42. 



elastic modulus at infinite frequency G'^. As both tran- 
sitions are approached, it tends to a finite value without 
showing any divergence, also in agreement with MCT. 

To make direct contact with experiments, we calculate 
G(v) = iujfftcj) — G'(w) + iG"{<J), where fj(ui) is the 
Fourier transform of Ccc (t) . To minimize numeric noise, 
we fit C aa (t) using functional forms that capture the fast 
short-time behaviour|29j and a long-time stretched expo- 
nential decay. Fig.0]shows the resulting G(oj) for the at- 
tractive glass transition case. Similar qualitative features 
are observed for the repulsive case. An incipient plateau 
at low frequencies in G'(u>) and a secondary maximum in 
G"(u>), leaving a minimum at intermediate frequencies, 
appear on approaching <tP . Both features are consistent 
with MCT predictions pa . \F\ and observed in experi- 
ments 

d 111 El The short-time behavior of C aa {t) 
causes the sharp increase in G'(ui) at w ~ 10 1 . 

In summary we have shown that in short range at- 
tractive colloids -q diverges following a power law as the 
transition points are approached, with the same expo- 
nent as the time scale T q , but different from that of Dq. 
On approaching (f> G , the slow-decaying part of C aa (t) 
can be time rescaled onto a master function which has 
the same shape as ® q (t) 2 at a specific q. For the case 
of repulsive glass, the inverse q value corresponds to the 
nearest neighbor distance, while, in the case of attrac- 
tive glass, it corresponds to distances comparable to the 



(short) range of attraction. Finally, the frequency de- 
pendence the elastic moduli is consistent, for both tran- 
sitions, with predictions based on MCT and with recent 
experiments on colloidal gels and glasses. 
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